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A generic chaotic eigenfunction has a non-universal contribution consisting of scars of short periodic 
orbits. This contribution, which can not be explained in terms of random universal waves, survives 
the semiclassical limit (when h goes to zero). In this limit, the sum of scarred intensities is a simple 
function of 77 = (/-l)^ 1 ^ A?) 1/2 , with / the de grees of freedom, /it the topological entropy 

and {Xi} the set of positive Lyapunov exponents. Moreover, the fluctuations of this representation 
go to zero as 1 / 1 In Tj, | . For this reasson, we will be able to provide a detailed description of a generic 
chaotic eigenfunction in the semiclassical limit. 

PACS numbers: 05.45.+b, 03.65.Sq, 03.20.+i 



Berry JjJ and Voros M proposed a semiclassical de- 
scription of chaotic eigenfunctions by considering the sur- 
face of constant energy as the unique classical invariant 
able to support them. This would be the case if the time 
required for the definition of individual eigenfunctions 
was infinite. However, the required Heisemberg time Tjj 
is finite for finite values of h, even though it goes to in- 
finity in the semiclassical limit (when h~ goes to zero). 

Recently, we have derived a semiclassical theory which 
describes individual eigenfunctions in terms of short pe- 
riodic orbits (POs) The number of POs is such 
that the sum of their periods is around Tg. Using these 
invariants, we will give a description of chaotic eigenfunc- 
tions with fluctuations going to zero with Ti. 

It is worth to emphasize that the fluctuations of the 
Berry- Voros description go to infinity in the semiclassi- 
cal limit (SL). But this statement demands explanation 
because it is common to find in the literature expresions 
like " ...this description is supported by Shnirelman theo- 
rem. ". For example, the peaks of the Husimi distribution 
of a chaotic eigenfunction have very strong fluctuations 
with respect to their classical ergodic measure; by assum- 
ing a random wave hypothesis they go logarithmically to 
infinity as h — ► |(|. However, each of the peaks oc- 
cupy a volume (the semiclassical volume of a quantum 
state) going to zero with h. Then, a smoothing on a 
region independent of %, a classical smoothing, is suf- 
ficient for washing out fluctuations in the SL. In this 
context, Shinirelman theorem |tJ may be expresed as fol- 
lows: a classical smoothing is sufficient for the elimina- 
tion of fluctuations ( with respect to the classical ergodic 
measure) of generic chaotic eigenfunctions in the SL. 

We are going to show that the highest peaks of the 
Husimi arc not completely random. On the contrary, 
they live (within an error of order yh in each direc- 
tion) on short POs, and their phases (considering in this 
case the corresponding Bargmann function |l) are con- 
nected semiclassically. That is, high fluctuations corre- 
spond with the scars of short POs introduced by Heller 
H (although the meaning of short is other; see Eq. (Eh). 



In this respect, the following commonly used expression 
is definitely wrong: "...scarred eigenfunctions do not con- 
tradict Shnirelman theorem because they are of null mea- 
sure. ". Actually, scar phenomena are of a semiclassical 
nature (they manifest on scales going to zero with K). 
Then again, a classical smoothing washes out these very 
rich structures in the SL. 

Bogomolny || and Berry (To) incorporated the idea 
of scars into the description of chaotic eigenfunctions. 
These authors concluded that chaotic eigenfunctions con- 
sist of a dominant universal contribution, decorated by 
scars of short POs; but this result was obtained after 
making an average over an energy interval that also elim- 
inates fluctuations. In our study, we consider individual 
eigenfunctions, without either a classical smoothing nor 
an energy averaging. Then, we arrive to the conclusion 
that a generic chaotic eigenfunction is scarred by a set of 
short POs which characterizes the state, and the sum of 
their intensities survives the SL. This fact stresses that 
localization on short POs is the signature of chaos in 
this limit. Moreover, we will give the mean value and 
dispersion of the scarred intensities showing that rela- 
tive fluctuations go to zero as l/|ln/x|, and providing a 
detailed description of a generic chaotic eigenfunction. 

We have developed a systematic semiclassical construc- 
tion of wave functions living in the neighbourhood of 
unstable POs. Reference || provides the construction 
of resonances without transverse excitations and ref. Q] 
applies the recipe to the Bunimovich stadium billiard. 
In ref. [jfj we construct resonances of a PO with trans- 
verse excitations, all at the same Bohr-Sommerfeld (BS) 
quantized energy. Then, scar functions are the result of 
a minimization of the energy dispersion in this basis. Fi- 
nally, ref. jll] gives the details for the construction of 
scar functions in the stadium billiard. Scar functions are 
the objects on which we are going to focus our descrip- 
tion. We are thinking about structures living along the 
manifolds up to the first homoclinic point. Then, it is 
not necessary to include interference effects, and the de- 
scription remains essentially simple. 
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Let 7 be an unstable PO of the system with period 
T 7 and Lyapunov exponent A 7 per unit time (consider 
for the moment a conservative Hamiltonian system with 
two degrees of freedom). Let 7 be the corresponding 
scar function with BS energy E 1 . Finally, let ip^ be 
the set of normalized eigcnfunctions of the system with 
eigenenergies E^. In rcf. j|] we have shown that the 
set of intensities 1^ = | (</> 7 |</? M ) | 2 semiclassically satisfies, 
E M h = !. E M E^I^ = £J 7 , and 

a 7 = y^7^~^)% - M 7 F/2 [12]. (1) 

In the SL, the universal dispersion Y goes to zero as 
27r/ 1 In h\ (expresions of T for finite values of h are given 
in ref. H). Then, the life time t 7 = Ti/a 1 of </> 7 di- 
verges logarithmically in the same way as the Ehrcn- 
fest time does (but they are different times). This ex- 
tremely low decay Jl3| is governed by a Gaussian law 
(|(<?!> 7 (0)|</i> 7 (£))| 2 = e~* ' t t), and the smooth part of the 
intensities I M (the strength function) results a Gaussian 
function. Then, defining e = (E — E^)pe (with pe the 
energy density), the mean value of the intensities at e is 

7(e) =e" e2/2 " 2 /V2W, (2) 

where n = a 1 pE is the mean number of eigenenergies 
contained in one energy dispersion. 

Another consequence of the Gaussian decay is that Eq. 
(|lj) is also valid for systems with / degrees of freedom if 
we replace A 7 by (J2 A 2 ) 7 ^ 2 , with {Ai} 7 the set of / — 1 
positive Lyapunov exponents of 7 (for an exponential de- 
cay the change would be A 7 by the K-S entropy Aj) 7 ). 

Now, we are able to establish a criterium in order to 
decide when a P.O is short. We propose a relation of the 
form T 7 < /3t 7 , with (3 a constant to be determined from 
the theory of short POs. This theory says that in the SL 
an eigenfunction is defined by POs of period lower than 
T = h^ 1 ln(Tf//iy) ||, with hx the topological entropy. 
Then, in the SL To would be equal to f3 times r (the life 
time of a generic scar function), with 

r = 7 r- 1 (^A 2 )- 1 / 2 |ln? l |. (3) 

{Xi} S ys is the set of positive Lyapunov exponents of the 
system. Finally, we say that a PO 7 is short when 

T 7 /r 7 < \p2/K j] " short PO condition", (4) 

where the classical invariant 

n^^/2(f-l)h T \J2^)lil (5) 

will play a central role below. Evidently, any PO is short 
at sufficiently high energies. 

Of all the intensities of a given scar function we are 
actually interested in those with the highest values (they 



are certainly related to the phenomenon of localization 
on short POs). We are going to study the properties 
of these high intensities in terms of a statistical model. 
The main purpose is to estimate their mean values and 
dispersions, and to provide a range in the spectrum where 
they live. 

The fluctuations of an intensity / at e, around its mean 
value 7(e) (see Eq. (||)), will be described as usual by 
a chi-squared distribution with one degree of freedom for 
systems with time reversal symmetry (appendix D treats 
systems without time reversal). Then, the probability 
of finding an intensity lower than / (the accumulated 
probability) is given, for a = I /1(e) > 1, by F t (I) = 

l-y/2/na e~ a/2 (l-a _1 +3a _2 H ) . In order to obtain 

the distribution for an arbitrary intensity (independent 
of its position in the spectrum), we make an averaging; 
if N is the number of intensities contributing to the scar 
function, F(I) = (1/jV) $ N (? /2 F t {I)de. This expresion is 
computed in terms of Gaussian integrals after the follow- 
ing expansion: e~ a / 2 = e~~ v e~ ye ~ ' l 2n (1 — ye 4 /8n 4 + ••■)) 
where y = \J^j2 nl (note that e 2 /n 2 = 0(y~ x )). Then, 

n/)=i-^^fi-A + j5L + ...y (6) 

W N y \ 8y 128y 2 J 1 ' 

Now, if x\ is the greatest of the intensities, X2 the 
second one and so on, the probability density of Xj is 
given by Eq. Jl6|). A simple estimation of Xj is derived 
from 1 — F(xj) — j/N. With this in mind, we define the 
random variable Zj by the relation 

l-F( Xj ) = e-»(j/N). (7) 

Combining Eqs. (|[) and (|?]), xj is given by 

y } = y/^flnxj ~ a - ln(a + 9/8) + b + b 2 /2, (8) 

with a = Zj + ]n(y/2n/j) and b = hx{a + 287/128)/(a + 
17/8). Equation (§) works for a > 1. 

Equation (|l7| ) gives the mean value ~zj and dispersion 
o~ Zj of Zj. With these, the mean value of yj is obtained 
from Eq. (||) by setting a (= ~z] + ln(y/2n/j)) in place 
of a and adding to the rhs the term c = a 2 j /2(a + 9/8) 2 

(because ln(a + 9/8) ~ ln(a + 9/8) — c). On the other 
hand, the dispersions of yj and Zj are equal to the leading 
order. Then, we arrive to the first conclusion: relative 
fluctuations go to zero in the SL as follows, 

oxjxi - l/[V7Mv^n/j)] = 0(l/\\nh\). (9) 

The next question is to know where Xj can be found. 
The probability p(e) of finding Xj near e is propor- 
tional to the probability density dF t /dI at Xj) that is, 
p(e) cx e -(w+ 1 /2)^/2n 2 ( 1 _ y je 4 /n 4 + ...). Then, x 3 is 
restricted to a range AEj ~ a 1 j yjyj + 1/2 around E 1 
which in units of 2nh/Tj (the distance between consecu- 
tive BS quantized energies) goes to zero as | ln?i|~ 3 / 2 . In 
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conclusion: scars of 7 accumulate in the vicinity of E 1 
in the SL. 

In the following, we are going to discuse a fundamental 
question. How many universal waves participate in the 
description of a chaotic cigcnfunction?. For maps the 
question is simple; there are N u waves with N u equal 
to the area of the map divided by 2irh~. However, for 
Hamiltonian systems it is not clear in general which one 
is the right, if any, Poincare surface of section. 

At classical level, we can say that during an evolution 
for a time lower than the period T m i„ of the shortest 
PO, it is impossible to extract non- universal information 
of the system. Using this time, Berry Jll| distinguishes 
universal aspects of the spectral rigidity from the non- 
universal ones. However, T m i n is not very useful because 
it is simply an estimation. Evidently, we need a time rep- 
resenting a short universal evolution but such that its de- 
termination is based on a deep knowledge of the system. 
For these reasons we propose the inverse of the topolog- 
ical entropy for this time. Then, if an eigenfunction is 
defined after an evolution equal to Th and we can do the 
travel in steps no greater than , where each step de- 
fines a universal wave, the number of required universal 
waves is 



(10) 



N u = T H h T . 



There is a way of verifying the accuracy of Eq. ( |10[ ) . In 
billiards, Birkhoff coordinates define the right Poincare 
surface of section (because all classical or quantal infor- 
mation is contained on the boundary). The area of the 
section is 2h~k£, with C the lenght of the boundary and 
k the eigenwave number. Then, N u = hC/it, and assum- 
ing that also Eq. ( |l0| ) is right, we arrive to the following 
expresion for the topological entropy h! T per bounce 



h' T = CI /it A, 



(11) 



with A the area of the billiard and / the mean length per 
bounce. We have verified Eq. (|ll|) within an error of 
2% in the stadium billiard with radius unity and several 
different areas. 

So far, we have described scar functions in the basis 
of eigenfunctions. Now according to the theory of short 
POs ||, eigenfunctions would be described equivalently 
in terms of scar functions. That is, imagine the spec- 
trum of scar functions given by all BS energies of all 
short POs. This means for a given PO, all BS ener- 
gies in those energy regions where 7 is short (following 
the criterium of Eq. (^)). Then, an eigenfunction ip u 
with eigenenergy E a is represented in the basis of scar 
functions by the intensities 7 7 = |(^|^> 7 )| , which are 
concentrated around E u in the spectrum of BS energies. 
The dispersion 07 = C£2(E 7 — E a ) 2 1^) 1 / 2 depends on the 
Lyapunov exponents of the orbits with BS energies in the 
vicinity of E u . However, in the SL there is a uniformiza- 
tion and the energy dispersion results independent of the 



position in the spectrum and equal to %/t, with r given 
by Eq. (H) . The smooth part of the intensities J 7 is given 
by Eq. (0), with 



ap E = T h /2itt, 



(12) 



and e = (E — E u )pe [@. If X\ is the highest of the 
intensities, X2 the second one and so on, the mean value of 
Xj is given by Eq. (0) (and discusion thereafter), and the 
relative dispersion by Eq. (||). Of all short POs, those 
with possibility of having intensity Xj satisfy | E a — E 1 \ < 

v/VVi + I/ 2 ; 

The following question is to decide a criterium for scar- 
ring. We will say that ip a is scarred by 7 if J 7 is greater 
than the greatest of the intensities provided by a random 
model of universal waves. The highest universal intensity 
is, to the leading order, = 2hi(N u )/N u @- Tnen > 
using Eqs. (|), (jTo)) and @, the condition 27 > 
reduces to 

yj > r\ " scarring condition" . (13) 
The number n scar of POs satisfying Eq. ([L3J) is given by 



/ 81n(l + IJ&rj) 
(1 + 477/5) 



(14) 



with 5 = (9 + v73)/2. This formula (and the next too) 
interpolates the behaviors for large 77 (obtained from Eqs. 
@ and (0)) and 77 going to zero (see appendix C). Finally, 
the sum of scarred intensities results 



(15) 



We emphasize that the unexpected results of Eqs. ( |l4| ) 
and (|l5|) depend decisevely on the use of scar functions 
in order to measure localization on short POs. On the 
contrary, by using wave packets in the transverse direc- 
tion to the motion (the so called vacuum states in ref. 
H), there result n scar /n — and ^2 Xj =0 in the SL. In 
this respect, the departure of universal behavior found in 
ref. Ea] has a weight going to zero in the SL. 

We stress that all formulae of statistical nature ver- 
ify an impresive agreement with numerical simulations. 
Moreover, Eqs. (^) and (^) were extensively verified in 
the stadium billiard. Finally, we present an example 
where all the ideas developed in the article have been 
tested. Figure 1) shows the decomposition of a very high 
excited chaotic state (plotted in configuration space in 
ref. ||). We found 9 scarred short POs (n scar ~ 8.5 
from Eq. (14)) providing a contribution of the 38% (Eq. 
( |l5| ) predicts 33%). Moreover, scarred intensities are in 
good agreement with predictions (see Fig. Id)). 

In conclusion, localization on short POs survives the 
SL and depends exclusively on 7/. Evidently, this local- 
ization will be strong in systems with few degrees of free- 
dom. In particular for / = 2, and assuming that A ~ hx, 
this localization results a universal property. 
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Appendix A. Let I\, I2, . ■ ■ , In be a set of indepen- 
dent random variables, with common probability density 
/(/), living in the range [0, 00). Let X\ be the great- 
est of the variables, xi the second one an so on. The 
joint probability density is given by p{x±, X2, ■ ■ ■ , xn) — 
N\f{xi)f(x 2 ) ■ ■ • f(x N ) ,{orxi>x 2 >...> x N ] other- 
wise p = 0. Then, the probability density of Xj results 
(after integration over the others variables) 



p(Xj) 



N\f( Xj ) F{ Xj ) 



N- 



(N-j)i(j-iy. 



(16) 



with F(x) = J Q X f(x')dx' the accumulated probability. 

Appendix B. We will derive the mean value and disper- 
sion of the random variable Zj defined in Eq. (^) . Assum- 
ing that zj = O(l), there results from Eq. (^) 1— F(xj) — 
£>(j/N), and then F(xj) N ~ j = exp [(N - j) lnF(xj)] = 
cxp [-AT(1 - F(xj)) + 0{j 2 /N)}. Moreover, AM/ (AT - 
j)\ = A J [1 + 0(j 2 /N)]. Using these approximations in 
Eq. (|l6|), the probability density of Zj is given by 

p( Zj )~fe-^+ e -^/(j-l)l. 

For instance, a numerical computation gives zj ~ 0.577 
and a Zl ~ 1.28. On the other hand, for large values of j 
we change to a new variable w = \fjzj, and expanding 
in powers of 1/y^J there results 
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with a — w 3 /6, b — w 4 (w 2 — 3)/72 and c = w 5 (5w 4 
45w 2 + 54) /6480. Finally, using this density we have 
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in excellent agreement with numerical data. 
Appendix C. For r\ going to zero, n s 



satisfies 



I(n scar /2) = /(«) 



'2/nri/n. Then, using Eq. 



there results n scar /n ~ 2 3 / 2 ^/ln(l/7y). 

Appendix D. The corresponding formulae for systems 
without time reversal are the following: yj = y/2nnxj ~ 
a - ln(a + 3/4)/2 + b + b 2 , in place of Eq. (§), with 
a = zj + lnU/2nn/j) and b = In (a + 31/16)/(4a + 5). 
Equation ( |i~3| ) is the same. Equation ( |l4| ) is replaced by 



r /n ~ e - " V27rln(l + l/77), 



and Eq. (15) by 



E 

3=1 



(i + v) 

a/3/4 + 77 



-2// 




FIG. 1. Linear Husimi density plots of: a) the state number 
141,755 of the desymmetrized stadium billiard with radius 1 
and area l + 7r/4, b) its non-universal contribution consisting 
of 9 scar functions, and c) its universal contribution consisting 
of 819 plane waves, d) The set of scarred intensities (dots) 
and the theoretical estimation curves x]±cr Xi ; horizontal line 
displays the value of the highest universal intensity. 
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